Phase Diagram of Cold Polarized Fermi Gas in Two Dimensions 
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The superfluid phase diagrams of a two-dimensional cold polarized Fermi gas in the BCS-BEC 
crossover are systematically and analytically investigated. In the BCS-Leggett mean field theory, 
the transition from unpolarized superfluid phase to normal phase is always of first order. For 
a homogeneous system, the two critical Zeeman fields and the critical population imbalance are 
analytically determined in the whole coupling parameter region, and the superfluid-normal mixed 
phase is shown to be the ground state between the two critical fields. The density profile in the 
presence of a harmonic trap calculated in the local density approximation exhibits a shell structure, 
a superfluid core at the center and a normal shell outside. For weak interaction, the normal shell 
contains a partially polarized cloud with constant density difference surrounded by a fully polarized 
state. For strong interaction, the normal shell is totally in fully polarized state with a density profile 
depending only on the global population imbalance. The di-fermion bound states can survive in the 
whole highly imbalanced normal phase. 

PACS numbers: 03.75.Ss, 05.30.Fk, 74.20.Fg, 34.90. +q 



I. INTRODUCTION 

The effect of Zeeman energy splitting h induced by a 
strong magnetic field between spin up and down elec- 
trons on Bardeen-Cooper-Shriffer(BCS) superconductiv- 
ity, which has been investigated many years ago [l], 
0, 0], promoted new interest in recent years due to 
the progress in the experiments of ultracold Fermi 
gases 0, 1 1 H 1|. The well-known result for 
weak-coupling s-wave superconductivity is that, at a 
critical Zeeman field or the so-called Chandrasekhar- 
Clogston(CC) limit h c = 0.707Ao where Ao is the zero 
temperature gap, the Cooper pairs are destroyed and a 
first order quantum phase transition from the gapped 
BCS state to the normal state occurs Further stud- 
ies showed that the inhomogcncous Fulde-Ferrell-Larkin- 
Ovchinnikov(FFLO) stateQ where Cooper pairs have 
nonzero momentum can survive above the CC limit up 
to /ifflo = 0.754Ao. However, since the thermodynamic 
critical field is much smaller than the CC limit due to 
strong orbit effect [l|, it is hard to observe the CC limit 
and the FFLO state in ordinary superconductors. 

Recent experiments on ultracold Fermi gas trapped in 
an external harmonic potential, serve as an alternative 
way to study the pure Zeeman effect on Fermi superflu- 
idity [1, 0, 0, 0, 0, Q. The atom numbers of the two 
lowest hyperfine states of 6 Li atom, iV-p and N± , are ad- 
justed to create a population imbalance or polarization 
P = (JVf - A r i )/(iV T + Ni), which simulates the Zee- 
man field h in a superconductor. The s-wave attrac- 
tion between the two hyperfine states is tuned around 
the Feshbach resonance to realize a strongly rather than 
weakly interacting Fermi gas. In three-dimensional case, 
the density profiles observed in experiments exhibit an 
unpolarized superfluid core in the center of the trap and 
a polarized normal gas shell outside 0,0)013) which jus- 
tifies that the ground state around unitary is a phase 
separation state(although there may be some dispute 
on the shell structure), predicted by early theoretical 



works [13 mUlEEE!. 

One of the theoretical interests in the study of po- 
larized Fermi superfluidity is to determine its phase 
structure in the whole interaction strength region, 
namely in the BCS-BEC (Bose-Einstein Condensation) 
crossover [l3|, [TJ]. The complete mean field phase dia- 
gram in coupling-imbalance plane as well as the critical 
Zeeman field h c and critical polarization P c are theoreti- 
cally predicted in three dimensional case (3D) [l3| . Since 
the s-wave mean field equations can not be solved analyt- 
ically in the whole coupling region, it is hard to determine 
a precise phase diagram in 3D even in mean field approx- 
imation. Recently, the quantitatively correct phase dia- 
gram in 3D has been obtained in quantum Monte Carlo 
calculations [Hi, 16]. However, the theoretical prediction 
of the phase diagram for a homogeneous system can not 
be directly examined in ultracold Fermi gas experiments, 
due to the effect of the external harmonic trap. To have a 
comparison with the experimental data, one should inves- 
tigate the phase diagram and the density profile in the 
presence of an external trap potential, using the same 
equation of state. In the case of 3D, the density profile 
can only be treated numerically [13] ■ 

While we have well understood the 3D phase diagram, 
the phase structure of polarized Fermi gases in low di- 
mensions promoted recently experimental and theoreti- 
cal interests. In one dimensional case the phase diagram 
is determined via exact solvable models [181 ] . In the two- 
dimensional case, while exact solvable models are lacked, 
the s-wave mean field equations can be solved analyti cally 
in the whole coupling parameter region 0, [2(| [U [22| . 
and the Fermi surface topology and stability condition 
are not trivial and different from those in 3D (23| . In 
this paper, we will determine the phase diagrams of a 
polarized Fermi gas in two dimensions, for both homo- 
geneous and trapped systems, and calculate the density 
profile of a trapped imbalance Fermi gas. Our results are 
totally analytical in the whole coupling parameter region, 
including the phase diagrams and the density profile. In 
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the final part of this paper, we also discuss the existence 
of di-fermion bound states in the polarized normal phase. 



II. BCS-BEC CROSSOVER IN TWO 
DIMENSIONS 

The BCS-BEC crossover problem in two dimensions 
has been widely discussed in the literatures p^l ii 
I22I ]. In this paper, we employ an effective 2D Hamilto- 
nian where the renormalized atom-atom interaction can 
be characterized by an effective binding energy (l9j. For 
a wide Feshbach resonance, the effective grand canonical 
Hamiltonian can be written as 



H = 



E 

<T=T,i 



d 2 r^(r) 



2M 



(J,-(T z h) ip a (r) 



U / d 2 r^(r)v4(r)^(r)V T (r), 



(1) 



where M is the fermion mass, fi is the chemical potential, 
U > is the contact attractive interaction, and a z = ±1 
correspond to a =|, j. We choose the unit h = 1 through 
the paper. The Zeeman field h can be created by either 
an external fieldp], [2, 3] or a population imbalance [3, Q. 
In the former case, the total particle number N is con- 
served, but the particles iVj and in the states f and 
I can transfer to each other[24j|, i.e., the chemical poten- 
tials for the two components are always the same, but 
the external field h induces an effective chemical poten- 
tial difference. In the latter case, and N± are both 
conserved, and the two chemical potentials can be ex- 
pressed as fj>i = [i + h and (i[ = fi — h. 

At finite temperature in 2D, the long range order is ab- 
sent and no phase transition can happen. At zero temper- 
ature, however, there do exist long range order [22j and 
one can safely consider phase transitions among different 
states. In this paper, we will study the phase diagrams 
at zero temperature in the BCS-Leggett mean field ap- 
proximation which is accepted to adequately describe the 
BCS-BEC crossover at T = [H|. 

In the balanced case with h = 0, the thermodynamic 
potential density of a uniform Fermi gas can be evaluated 
as ['221 



A 2 
U 



d 2 k 



(2) 



with the definition of particle energies = v/^ + A 2 
and £k = k 2 / (2M)—/i and the superfluid order parameter 
A = —U{ij},ij}^). In the dilute limit, the UV divergence 
in the expression of f2 can be eliminated via introducing 
the two body scattering length. 

While the two-body bound state in 3D forms only at 
sufficiently strong attraction where the s-wave scattering 
length diverges and changes sign, the bound state in 2D 
can form at any arbitrarily small attraction [26|. For an 
inter-atomic potential described by a 2D circularly sym- 
metric well of radius Rq and depth Vq, the bound-state 



energy e B is given by e B ~ l/(2Mi? 2 ) exp [-2/ (MV0R 2 .)} 
with VoRq — » 0. As a consequence, the solution of the 
BCS-BEC problem in 2D is much simpler than that in 
the case of 3D in terms of special functions [20]. It is 
shown that the existence of the two-body bound state in 
vacuum is a necessary (and sufficient) condition for the 
Cooper instability [l9j. To regulate the UV divergence in 
f2, we introduce a high energy cutoff A = k 2 v /(2M) in 
the integral. The momentum cutoff kA corresponds to 
the inverse of the range r$ of the interaction potential. 
Due to the energy independence of the density of states 
in 2D, after performing the integration over k one obtains 



ft = — 



A 2 

U 



MA 2 

47T 



ln 



A - n + y/(A - + A 2 



vV + A 2 



A - /i 



/' 



A - n + y/(A - ^) 2 + A 2 V^ 2 + A 2 - /i 



(3) 



In this paper we consider a dilute Fermi gas with effec- 
tive interaction range ro — > 0. Taking large enough cutoff 
A and small enough attraction U , we can introduce a 2D 
two-body binding energy [22] to replace the cutoff in this 
limit. 



, 4tt \ 

£r = 2 A exp , 

1 1 MU J ' 



(4) 



which does not include any many-particle effect. With 
the binding energy, the cutoff dependence can be elimi- 
nated in the dilute limit with A — > 00 and U — > but 
finite £b- We obtain in this limit 



ft= MA /yZ±A^ 



1 



4tt 



vV + A 2 - n 2 



. (5) 

The above procedure is equivalent to directly substitut- 
ing the coupling constant U by the 2D bound state equa- 
tion 



1 r d 2 k 



U J (2tt) 2 k 2 /M + 6 B ' 



(6) 



The BCS-BEC crossover phenomenon in 2D can be 
observed by solving the coupled gap and number equa- 
tions, namely dfl/dA = and n = —dfl/dii. Defining 
the Fermi energy ef = ttu/M in 2D, the gap and number 
equation can be analytically expressed as 



vV + A 2 - M = e B , ^ 2 + A 2 +M = 2e F , (7) 
respectively. Their solution takes a very simple form[l9| 



A = V2£b£f , Mo = £f-y, (8) 

or rewrite it in terms of a dimensionless quantity rj = 
£b/£f, 



An 



27], 



e F 2 



(9) 
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One sees very clear that the chemical potential decreases 
with increasing coupling or decreasing density, which in- 
dicates a BCS-BEC crossover. The Chemical potential 
changes sign at ea = 2ep. To understand the physical 
significance of these simple results, we consider two lim- 
its. For very weak attraction (or high density), the two- 
particle binding energy is extremely small, i.e. ee < £f, 
and we recover the well-known BCS result with strongly 
overlapping Cooper pairs. In this limit we have the chem- 
ical potential /Lto — £f an d the gap function A -C e-p. For 
the opposite limit of very strong attraction (or low parti- 
cle density), we have a deep two-body bound state with 
ee > £f, and the system is in the BEC region with com- 
posite bosons. In this limit the chemical potential takes 
fio ~ — cb/2. It should be kept in mind that in the local 
pair regime (/iq < 0) the fermion excitation gap -E gap in 
the quasi-particle excitation spectrum is not Ao (as in 
the case /j,q > 0) but rather y/ ^ + Ag. 

In ultracold Fermi gas experiments, a quasi-2D Fermi 
gas can be realized by arranging a one-dimensional opti- 
cal lattice along the axial (z) direction and a weak har- 
monic trapping potential in the radial (x-y) plane, such 
that fermions are strongly confined along the z direc- 
tion and form a series of pancake-shaped clouds [2?], HI] . 
Each such cloud can be considered as a quasi-2D Fermi 
gas when the axial confinement is strong enough to turn 
off inter-cloud tunnelling. The strong anisotropy of the 
trapping potentials, namely u) z 3> lu where uj z (lo) is 
the axial(radial) frequency, allows us to use an effec- 
tive 2D Hamiltonian to deal with the radial degrees 
of freedom] 28]. The effective binding energy ee = 



ftw z exp [47ra^/t7p ff (a s , a z )] is related to the energy scale dA 



tko z and the 3D s-wave scattering length a s , where a z is 
defined as a z = y/h/ (mui z ) and the quantity /7° ff (a s , a z ) 
defined in [28j carries the dependence on 3D scattering 
length. By adjusting the 3D scattering length a s and/or 
the axial frequency u> z , a quasi-2D BCS-BEC crossover 
can be realized. 



consistently from the gap equation 

1 Q(E Z - h) 



A 



dz 



2z + e B 



2E* 



(11) 



The step function 0(x) in this paper is defined as 0(x) = 
for x < and Q(x) — 1 for x > 0. 

Unlike the 3D casel3|, the EOS in 2D can be an- 
alytically obtained. At fixed /i and h, we have three 
possible phases: the unpolarized superfluid phase (SF), 
the polarized normal phase (N) and the polarized su- 
perfluid phase or Sarma phase (S) 0. The phase SF 
corresponds to the solution A(/i) = ^/cb^b + 2//) in the 



region h < E g = y/ /i 2 0(— /i) 
can be evaluated as 

VsFiv) =c(fJ, 



A 2 (/i) , and the pressure 



£b 
2 



(12) 



which does not depend explicitly on h. The total number 
density n — + n± and the magnetization m = w\ — 
can be expressed as 



«sf(a«) 



2c i 



(m+y), m SF (/i)=0. (13) 



The polarized superfluid phase or Sarma phase (S) corre- 
sponds to the solution in the region h > E g . This phase 
can be ruled out from the positive secondary deriva- 
tive m 



d 2 V 



fe(l + e(M)) _ A*6(M) 
Vh 2 - A 2 yfv 2 + A 2 



- 1 



>0, (14) 



which means that the Sarma phase is always unstable for 
any coupling in 2D. 

The polarized normal phase corresponds to the so- 
lution A = 0. The pressure takes the form of non- 
interacting Fermi gas, 



III. EQUATION OF STATE 

We now turn on the Zeeman splitting h ^ 0. To de- 
termine the superfluid phase diagrams and calculate the 
density profile for imbalanced Fermi gas in a harmonic 
trap, we first establish the equations of state (EOS) for 
various phases in grand canonical ensemble [13j. In the 
BCS-Leggett mean field theory, the pressure V = —Q as 
a function of fi and h can be evaluated as[l3l] 



dz 



E z 



z + fi - 



2z + e B 



{E z -h)9(h-E z 



(10) 



with E z = y/(z - n) 2 + A 2 and c = M/(2n). We have 
set h > without loss of generality. The superfluid or- 
der parameter A(/j,,h) = —U(ip^ip-^) is determined self- 



TMm, h) = - [(n - h) 2 G(fi -h) + (fx + h) 2 B((i + h)] , 

(15) 

where the case /j, + h < corresponds to the vacuum 
without atoms. For fj, + h > 0, the total number density 
and the magnetization read 

71n(a«, h) = 2c/i© (fi — h) + c(fi + h)Q (h — /i) , 
m N (p,h) = 2chB {fi - h) + c(p + h)G (h - fj,) . (16) 

The cases fi > h and fi < h correspond to the par- 
tially polarized (Npp) and fully polarized (Npp) normal 
phases respectively. Since we treat the superfluid and 
normal phase in mean field approximation, the normal 
phase is considered as a non-interacting gas. In fully po- 
larized case, this is correct since only s-wave interaction 
is considered. However, in partially polarized case, the 
interaction may be important in some coupling parame- 
ter region, like the finding around the unitary region in 
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3D [l6|, [2^ . Including fluctuations, which can not be 
treated analytically even in 2D, is necessary for a more 
realistic study. 

Since the polarized superfluid phase is always located 
at the maximum of the thermodynamic potential, there 
exists at fixed /1 a first order quantum phase transition 
from the SF phase to the normal phase when the Zee- 
man field h increases. The critical value h c is determined 
by the condition PsfO-0 = ^n(a*j h c ). The analytical 
expression for h c can be written as 



<"B 



)e( M -M 
e(h -fi). 



(17) 



Equivalently, for a given h, SF-N phase transition hap- 
pens when the chemical potential /x becomes less than 
the critical value 



fj, c (h) 



e B 4 
V2h-e B 



V2 



e(h-ho) 

@{ho-h), (18) 



where ho = (\/2+ l)e B /2 is determined by the equation 
^0 = Mc(/ l o)- We can easily prove that h > ho(h < ho) is 
equivalent to the condition /i c > h(fi c < h). 

The grand canonical phase diagram in the fi — h plane 
is shown in FigfTJ The analytical expressions for the 
phase boundaries can be obtained from the above ex- 
pression for h c . The SF phase, Nfp phase and the vac- 
uum meet at the point (n,h) = (— es/2, eb/2), while 
the three phases SF, Npp and Npp meet at (/j, h) = 
((V2+l)eB/2, (V2+l)eB/2). The grand canonical phase 
diagram is of great help for us to understand the density 
profile in a harmonic trap. 




FIG. 1: The grand canonical phase diagram in the fjb—h plane. 
(i and h are scaled by the binding energy ee- 



IV. HOMOGENEOUS FERMI GAS 

In this section we determine the phase diagram of 
the homogeneous system. Since the total atom number 
N = Nt + N\ or equivalently the total atom density n 
is fixed, the chemical potential /1 is not a free variable in 
the canonical ensemble and should be determined by the 
number conservation. One may distinguish two different 
cases: (l)The Zeeman field h can be experimentally ad- 
justed by using Raman detuning [24]]; (2)The atom num- 
ber for each species, iVj and N±, can be adjusted[l, [f|. 
Since the phase structure should be essentially indepen- 
dent of the ensemble we choose, we firstly discuss the 
phase diagram using h as tunable parameter, and then 
translate it into the case where the global polarization 
P = (TV-f- — A r i)/(A^T + N i) is directly adjusted. 



A. Critical Zeeman Fields 

We now consider the problem: When does the super- 
fluidity disappear when a Zeeman splitting h is turned 
on? The total density n and the Zeeman field h are ther- 
modynamic variables, and the free energy of the system 
should be defined as J-{n, h) — fxn—V. Since n = Mcf/tt 
is fixed, we will write 3-{n, h) = J-(h). At nonzero Zee- 
man field h, the solutions of the coupled gap and number 
equations corresponding to the above three homogeneous 
bulk phases can be analytically solved: 

I. AsfCO = A and fJ.sp{h) = Ho hi the phase SF. The 
solution exists in the region < h < Ao for rj < 2 or 
< h < e F + e B /2 for r\ > 2. 

II. A N (/i) = and fi N (h) = e F Q{e F -h) + (2e F -h)e(h- 
e F ) in the phase N. The first and second term corre- 
spond, respectively, to partially and fully polarized nor- 
mal phase. 

III. A s (h) = ^A (2/i- A ) and fi S (h) = e F - 
A|(/i)/(4ep) in the region of A /2 < h < A and 77 < 2 

and Ag(/i) = y e B (2/i — ee) and Hs(h) = 2e F — h in the 
region of 6b/2 < h < e F + e B /2 and 77 > 2 in the phase S. 
There are two gapless Fermi surfaces at r\ < 2 and only 
one gapless Fermi surface at 77 > 2. 

The polarized superfluid phase or the Sarma phase, 
which is a gapless superfluid, is again an unstable state at 
any coupling, directly from the reentrance phenomenon 
(three solutions of A at fixed h) , in contrast to the case 
in 3D where it becomes the stable ground state in the 
strong coupling BEC region [H,[i3]- This is an important 
difference of the Fermi surface topology and the stability 
condition between 3D and 2D cases [23|]. Explicitly, the 
free energy (density) T{h) = fx(h)n — V(fi(h), h) in the 
three homogeneous bulk phases reads 



r SF {h) 



c(4 

c[(e| - ti 2 )Q{e F - h) 
2(el-e F h)G(h-e F )], 
c[2(el - e F h) + h 2 + (h 



e B /2) 2 ]e(?7-2) 
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c[(4 



e F e B ) + (A -/ l ) 2 ]e(2-ry). (19) 



It is easy to see that the polarized superfluid phase has 
always higher free energy. If there exist no other possible 
phases, a first order quantum phase transition from the 
phase SF to the phase N will occur at a critical Zecman 
field h c determined by Tsw{h c ) = T-^{h c ). We find h c — 
,/fjeF = A0/V2 for 77 < 1 and h c = |(1 + rj)ep for 77 > 1. 
It is interesting to note that the relation h c = Aq / \/2 at 
77 < 1 is only an approximate result at weak coupling in 
3D 0,0. 

If there are only the two bulk phases SF and N, we have 
only one CC limit at which the first order phase transi- 
tion occurs, and the experimentally observed phase sep- 
aration (PS) will be hidden in the 77-/1 phase diagram. 
However, since the total atom density n is fixed, un- 
like the grand canonical ensemble, we should consider 
possible mixed phases constructed via the Gibbs phase 
equilibrium condition. Here we will neglect the interfa- 
cial energy [3(j> since for a macroscopic phase separation, 
this energy contribution is subdominant in the thermo- 
dynamic limit. From equation (|14[) . the only possibility 
is the SF-N mixed phase. When the phase separation is 
favored in a region h c \ < h < h C 2 , the chemical potential 
fips is different from /usf and /un, it should be determined 
by the phase equilibrium condition "Psf(/j) = V^^fijh), 
which leads to 



Mps(^) 



_ £2. 

e B 4 

V2h-e B 
2-V2 



®{h-h ) 

e {h - h) , 



(20) 



where ho — (\/2 + l)es/2 satisfies the equation ho = 
HPs(ho). Since the chemical potential fips is determined 
by the condition 'Psf(m) = P^{^,h), it is equivalent to 
the critical chemical potential fi c (h) in the grand canon- 
ical ensemble. The cases h > ho and h < ho indicate, re- 
spectively, the mixed phases with partially polarized nor- 
mal bubbles (SF-Npp) and fully polarized normal bub- 
bles (SF-Npp). The volume fractions of the phases SF 
and N in the phase separation, denoted by x and 1 — x 
respectively, are determined by the number conservation, 
n = x(h)nsF(fJ,ps,h) + [1 — xQi)]njq(fips, h). Using the 
expressions fllUl) . (fTO|) and (|2"P|) for /Jps,«N and nsF, we 
find 



x(h) 




h 2 " 

— )Q(h- h ) 

e B 



2h - e B 



V2-l]e(ho-h). (21) 



We now determine the region of the mixed phase, i.e., 
the lower and upper critical fields h c \ and h C 2 [l3j]. In 
the grand canonical ensemble with fixed chemical poten- 
tial n, we have only one critical field h c (fi) determined 
by the condition Vsp(fJ-,h) = V^{^h), and the signal 
of SF-N phase separation is denoted by the first order 



phase transition line in the fi — h phase diagram. In the 
standard BCS-BEC crossover problem, the total atom 
number N rather than the chemical potential /i is fixed, 
and the CC limit splits into two values h c \ = h c (fisp) 
and h C 2 = h c (n^) [131 ]. The mixed phase links continu- 
ously the phases SF and N with /ips — (J-sf at h = h c \ 
and /ips = 1-J.n at h = h C 2 and ensures < x < 1 with 
x(h c i) = 1 and x(h C 2) = 0. The critical fields h c \ and 
h C 2 are explicitly given by 



hd 



= ^F^Jr](l- 6(771 - 77), 

+ e F (72- 1 + 6(77- 771), 

- ^(1 + 1)6(772-77), 

+ e F (2 -V2 +|)e(77 -772), 



(22) 



where 771 = 2 - y/2 ~ 0.586 and 773 = 2(y/2 - 1) ~ 
0.828 are determined by /is F (ft.o) = ^0 and ^(ho) = ho, 
respectively. There is always the relation h c \ < h c < 
h C 2, and the splitting disappears in the weak coupling 
limit ?7_— > which recovers the well known result shown 
in [H, On the other hand, the splitting keeps as a 
constant (3 — 2y/2)ep ~ 0.172ep at strong coupling 77 > 

The final step is to prove that the SF-N mixed phase 
has the lowest free energy in the region h c \ < h < h c2 . 
Using the analytical expressions for x(h) and fips(h) 
as well as the EOS for the phases SF and N, we can 
evaluate the free energy in the mixed phase defined by 
fps(h) = Hpsn-x(h)VsF([i>PS,h)~[l-x(h)]'P N (iips, h). 
The difference between .Fps and Tsf and between ^ps 
and can be explicitly expressed as 

^ps(fc) ~ Fs?{h) = -t B 2 c(h 2 - h 2 cl ) 2 e(h - h ) 
-{V2 + l) 2 c(h-h cl f@(ho-h), 



Fps(h) - Mh) 



c{h z - hi 2 ) 2 e(h - h ) 



-(V2 + l) 2 c(h - h c2 ) 2 Q(h - h).(23) 

The above expressions show explicitly that the mixed 
phase has really the lowest free energy in the region h c ± < 
h < h C 2- While in 3D the conclusion that the mixed 
phase corresponds to the lowest free energy is analytically 
proven in the weak coupling limit [10j , our result here in 
2D is for any coupling. 

The SF-N mixed phase has a nonzero global polariza- 
tion P since the normal bubble is polarized. From the 
definition P = (7V T - N l )/(N- t + AT;) we find 



P(h) = [l-x(h)] 



m N (/Jps, h) 



(24) 



Using the expression for x(h) and /ips, we have 



P(h) = 2 ^ 2 -^W- M 



£ f£ b 



{V2 + l) 2 (h-h cl ) 



CF 



Q(ho-h). (25) 
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The global polarization is zero at h 
creases with h. 



hri and then in- 



B. Critical Polarization 

Finally, we convert the above result into the one where 
both JVf and N± are fixed and the exchange between par- 
ticles in states | and | is forbidden, corresponding to 
recent experiments on ultracold Fermi gas with popu- 
lation imbalance 0, The free energy density in this 
case should be defined as !F[n-\, n{) = /i| n T + — P 
or Tin, m) = fin + mh — V. The possible phases with 
nonzero global polarization P are the normal, Sarma and 
SF-N mixed phases. Since the phase structure should be 
essentially independent of the ensemble we choose [13|, 
we do not need to compare again the free energies of the 
three phases [Toj. Since P(h c i) = and P(h) increases 
with h, we conclude that the ground state is the unpolar- 
ized superfluid at P = 0, and SF-N phase separation be- 
comes energetically favored for < P < P c . The critical 
polarization P c where the superfluid bubble disappears 
completely is the global polarization at h C 2, 

P(h c2 ) = —Q(m -v) + ®{v - m) 



Pr 



r?(i + j)®(v2 



??) + Q(r) - r) 2 ). 



(26) 



The critical polarization increases from P c = at r; = 
to P c — 1 at 77 — r\2 an d then keeps as a constant P c = 1 
at strong coupling 77 > 7/2 



C. Phase Diagrams 

FigfJ] summarizes the above analytical results. The 
phase diagram in the rj — h plane is shown in the upper 
panel. The partially and fully polarized normal phases 
TVpp and -/Vpp are separated by the dashed line h/ep — 1 
which ends at 77 = 772. The two solid lines indicate the 
lower and upper critical Zeeman fields h c \ and h C 2 with 
the two phase separations PS-I and PS-II in between. PS- 
1 (SF-Npp)and PS-II (SF-N F p) are the mixed phases of 
superfluid and normal gas with JVpp and iVpp, and they 
are separated by the dotted line h/e-p = (v2 + 1)77/2 
starting at 1} = r/i and ending at 77 = 772. The phase 
diagram in the 77 — h plane can be easily converted into 
the one in the 77 — P plane shown in the lower panel, by 
taking the fact P(/i cl ) = and P(h C 2) = P c - The critical 
polarization P c — \J 77(1 + 77/4) (solid line) increases from 
P c = at 77 = to P c = 1 at 77 = 772 and then keeps 
as a constant P c = 1 for 77 > 772. The phases SF and 
iVpp are now located at P = and P = 1 respectively. 
The dotted line which separates PS-I from PS-II becomes 
P = (4 + 3V2)?7/2 - (V2 + 1) in the 77 - P plane. 

The above analytical results show that, to correctly 
calculate the critical polarization P c and the phase dia- 
grams, one should treat the mixed phase carefully [l3l |. 
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FIG. 2: The phase diagrams in the planes 77 — h (upper panel) 
and n — P (lower panel), h is scaled by the Fermi energy. 
SF means unpolarized superfluid, iVpp and N-pp indicate the 
partially and fully polarized normal phases, and SF-Npp and 
SF-Nfp are the mixed phases of superfluid and normal gas 
with iVpp and Npp. 



Some other methods taken in literatures may lead to 
quantitatively incorrect results. For instance, the method 
of stability analysis will result in an incorrect critical po- 
larization (see also the comments in [Hj]). With this 
method, one first solve the mean field gap and number 
equations for the Sarma phase and then analyze the sta- 
bility of this phase. If it is ap plie d to the 2D system, the 
critical polarization becomes [3 ll ] 



Pr. = 



2e F 



1 
2" 



(27) 



which is the maximum polarization of the unstable Sarma 
phase and deviates significantly from our result P c = 
77(1 + 77/4). Especially, our critical polarization grows 
up to unity at 77 ~ 0.828, but the result (|2T|) becomes 
unity at 77 = 2. On the other hand, if one takes only 
the phases SF and N into account but neglect the phase 
separation, there will be only one critical field h c where 
the polarization jumps from to h c / e-p = <Jfj. 
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V. BOUND STATE IN POLARIZED NORMAL 
PHASE 

In recent experiment on highly polarized normal phase 
in 3D unitary Fermi gas[8j, it is found that while the su- 
perfluidity disappears completely, full pairing of minor- 
ity atoms always exists, which indicates that the fermion 
pairing may be easy to occur in the presence of polariza- 
tion. It is well known that the bound state in 2D can 
form at arbitrary small attractive interaction [26[, which 
is quite different to the 3D case. It is natural to ask: Do 
the di-fermion bound states exist above the upper critical 
field h C 2 or critical polarization P c ? 

In this section, we study the spectrum of bound states 
in the highly polarized normal phase. In the Green func- 
tion method, the energy lu of the bound states with zero 
total momentum in this case is determined through the 
equation [22j 



1 



i - e(/i T - z) - e(/i| - z) 



2z + e B 



2z — 2/i — u) 



dz ^ i ^± '- =0. 

(28) 

In the vacuum with /_t = h = 0, it self-consistently gives 
the solution lu = — ee- In general case with medium ef- 
fect, the bound states can survive when the above equa- 
tion has real solution of lu. 

The integration in the above equation can be analyti- 
cally worked out, and finally we obtain 



, lu + 2u lu + 2h 

In - + 6(/i-ft.)ln- 



2/i 



-e( M +/ l )ln^-|^ = 0. 

LU + Zfl 



(29) 

In the partially polarized normal phase, we have fi = ep, 
the spectrum equation becomes 



lu a + e B LU + 2e F e B - 4/r = 

which has real solutions 
1 



LU 



for 



J(h) = e| + I6h 



i6 F > 0. 



(30) 



(31) 



(32) 



In the fully polarized normal phase, one finds that the 
spectrum equation directly gives a real solution lu — 
2h — 6b- However, this solution is unphysical since we 
always have lu > 0. This can be well understood when 
we consider the fact that there exist only f particles in 
this phase. 

In the balanced normal phase with h — (note that 
the true ground state in this case is the superfluid phase) , 
the bound states remain stable only at strong enough 
coupling r] > 8 or equivalently low enough density e-p < 
cb /8 [22j , which indicates that the Fermi sea or medium 
effect disfavors the formation of bound states. One may 
simply think that, the presence of a Zeeman splitting 
will further destroy the bound states. However, this is 




r|=e B /e F 



FIG. 3: The excitation gaps for the di-fermions and di-holes 
as a function of rj at fixed polarization P = 0.8. The curves 
are meaningful only for r\ < 0.56, since the ground state is 
not a normal phase for r\ > 0.56. 



not true. Since the condition J(h) > is easier to be 
satisfied at h ^ 0, the bound states in highly polar- 
ized normal phase are easier to survive than in the bal- 
anced Fermi sea. In the whole partially polarized normal 
phase which exists in the region < rj < rj2 in Figj^l we 
have 6F\^r](l — 77/4) < h < ep. Analyzing the condition 
J{h) > 0, we conclude the bound states can exist in the 
whole Npp phase in FigO Especially, they can survive 
at high polarization even in the weak coupling limit, as 
pointed out in [33| in the case of 3D. 

There exist two real solutions for lu in the Npp phase. 
The negative {lu < 0) and positive (lu > 0) solutions in 
(|3"Tj) correspond to the excitation gaps E^f and E^h for 
the di-fermions (df) and di-holes(dh) respectively [33), 



E di 



2 



(?/ + vV - 877 + 16P 2 ) , 



E dh = S (y r? 2 - 87/ + 16P 2 - 7/ 



(33) 



In Fig|21 we plot the excitation gaps for di-fermions and 
di-holes at a fixed polarization P — 0.8. We found that 
the symmetry in the spectrum(£ , df = ~Pdh) holds only 
at weak coupling. 



VI. DENSITY PROFILE IN A HARMONIC 
TRAP 

We have determined the phase diagram for homoge- 
neous system. However, the phase structure can not be 
directly examined in ultracold Fermi gas experiments, 
due to the effect of the external harmonic trap. To jus- 
tify the theoretical prediction for homogeneous system, 
one should calculate the corresponding phase diagram 
and the density profile in the presence of an external 
trap potential, using the same equation of state. In this 
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section, we will calculate analytically the density profile 
of an imbalance Fermi gas in a 2D isotropic harmonic 
trap potential V(r) = \Mio 2 r 2 . The frequency u> here is 
different from the energy of the bound state defined in 
Section V. 

The effect of a harmonic trap can be treated in the 
local density approximation(LDA). In the frame of LDA, 
the system is approximately taken to be uniform but with 
a local chemical potential given by 



M( r ) = Mo - l;Muj 2 r 2 , 



(34) 



where Mo is the chemical potential at the center of the 
trap and is the true chemical potential(a Lagrangian 
multiplier) still enforcing the total atom number N = 
iVf + N±. Since iVj and N± are both conserved, the 
spatially-varying spin-up and spin-down local chemical 
potentials can be expressed as Mt( r ) = M( r ) + h an d 
Ml( r ) = M( r ) ~ h hi terms of the averaged chemical po- 
tential fj,(r) and Zeeman field h. 

To calculate the density profile, namely the atom den- 
sity of the spin-up and spin-down states as a function 
of the radius r, n-f(r) and ni(r), or equivalently the to- 
tal density n(r) = H|(r) + ni(r) and the magnetization 
m(r) = ri]-(r) — ni(r), one should know the equation of 
state, n a (r) = n a (fi(r),h) with a =\,i- Using the EOS 
calculated in Section III, we can determine /iq and h from 
the known total particle number N and the global polar- 
ization P = (JV t - N l )/(N- l + N{), 

N = 2tt J rdrn(r), PN = 2n J rdrm(r). (35) 

Let us firstly consider a non-interacting system with 
balanced populations, N-\ — N^, which can help us to 
define the Fermi energy e F for trapped 2D system, e F is 
defined as the chemical potential Mo at the center of the 
trap for non-interacting gas. The density profile is also 
balanced, n^(r) — nj,(r), and is given by 



£F 



(36) 



which vanishes at the so called Thomas-Fermi radius 



i?T — 



2e F 
Mlo 2 



The total density N is then given by the integral 



V = 2k I rdr2c ^e F - 



(37) 



(38) 



We find that the Fermi energy in 2D is e F = yNhiu, in 
contrast to the result e F = {QNY^fouj in 3D (note that 
we have recovered Ti in these expressions). 

We then turn to an attractive Fermi gas. For balanced 
populations, the ground state is a superfluid state, and 
the density profile can be obtained from equation (Q2 

n l( r ) = n l( r ) = c (mo + y - ^Mu 2 r 2 



(39) 



Comparing with the non-interacting gas, we find no dif- 
ference between the normal and the superfluid states. 
The chemical potential at the center of the trap reads 



Mo 



e F 



£B 
2 ' 



(40) 



This relation is exactly the same as in the homogeneous 
case[lj|. The order parameter profile A(r) is given by 



A(r) 



1 - r 2 /i?^, Ao = V2^F- 



(41) 



For a system with population imbalance, iVj ^ 
we should have h ^ 0. By comparing Psf with Vn, 
a first order phase transition from the phase SF to the 
phase N occurs for a given /i when the Zeeman field h 
becomes larger than the critical value /i c (m)i or equiva- 
lently speaking, for a given h the SF-N phase transition 
happens when the chemical potential /i becomes less than 
the critical value n c (h). In LDA, the phase behavior as a 
function of chemical potential m is translated into a spa- 
tial cloud profile through mM- The critical phase bound- 
ary Mc corresponding to the critical radius r c is defined 
by 



M c = fi(r c ) = /i 



Mu 2 r 2 



(42) 



at which the states SF and N have the same pressure. 
Thus, at fixed h, any region of the system which satis- 
fies M( r ) > Mc is in the state SF, while a region which 
satisfies fi(r) < fi c will be in the state N. Since /i(r) de- 
creases with increasing r, it is clear that the high density 
superfluid region will be confined in the center of the 
trap, and the low density polarized state N is expelled 
to the outside. The shell structure with radius r c of the 
SF-N interface is a striking signature of phase separation 
in a trap. The superfluid core will disappear when the 
population imbalance P becomes larger than the critical 
value P c which is determined by the equation r c = + or 
Mo = Mc- 

We should have two types of shell structure corre- 
sponding to the cases h > ho and h < ho. For the case 
h > ho, we have n{r c ) = Mc = h 2 / e B — £b/4 > h, which 
means that there exists a shell of partially polarized nor- 
mal gas in the region r c < r < ro, with ro given by 
M( r o) = h. We call it the phase PS-I. Thus we have the 
following density profile 



n(r) 



2c (mo t- — 



$f - ±Muj 2 r 2 
2c (/io - I 
c (mo +h- 



^Muj 2 r 2 ) 
Mw 2 r 2 ^ 



<r <r c 
r c <r <r 
r <r < R 



and 



m(r) 



(43) 



(Mo 





2ch 
\Mu 2 r 2 ) 



< r <r c 
r c <r <r (44) 
r <r < R 



where R = ^/2(/i + h)/Mui 2 is the edge of the cloud. 
After some algebra according to the equation (|3"5"|) . Mo is 
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simply given by fio = cf — £b /2 as in the balanced case 
and h is solved from the cubic equation 



2h 



£b 
4 



rf Fi 



(45) 



where ep = ^/Nhio is the Fermi energy defined in (|38|) . 
From the condition h > ho which ensures rg > r c , we 



have Pep > 2/i 2 , which leads to the relation 



P>Pn = 



2^2 



-'7 



(46) 



The critical polarization P c is determined by the condi- 
tion hq = /i c . A simple algebra gives 



Pr 



(2 - 77)1/77 



T 

T 



< 77 < 771 



(47) 



with 771 = 2 - y/2 ~ 0.586. Note that both P and P c 
reach unity at 77 = 771 , they are the two boundaries of the 
phase PS-I in the 77 — P plane. 

For the case h < ho or P < Pq, we have /i(r c ) = /_t c = 
(V2/1 - e B )/(2 - \/2) < ft, which means that the normal 
gas shell outside the superfluid core is fully polarized. 
The density profile reads 



n(r) 



f 2c (mo + f - I 



Mw 2 r 2 



< r < r c 
r r < r < R 



and 



m(r) = 



< r < r c 



c(/i + ft-- |MwV) , r c <r<R 



(48) 



(49) 



After the integration in equation (|35[) . we still have 
ep — £b/2 and /i is explicitly given by 



h=(V2- l)VPe F 



(50) 



One can easily check that the condition h < ho is equiv- 
alent to P > Po, and we have P c = 1 for 77 > 771. 

Fig[4] summarizes the the coupling-imbalance phase 
diagram for two-dimensional imbalanced Fermi gas in 
a harmonic trap. The critical polarization P c = (2 — 
i])^/t] — 77 2 /4 (solid line) increases from P c = at 77 = 
to P c = 1 at 77 = 77! ~ 0.586 and then keeps as a constant 
P c = 1 for 77 > 771. The dashed line, analytically given 
by P = (3 + 2v / 2)?7 2 /2, separates the two types of phase 
separation, PS-I and PS-II with different shell structure. 
In the phase PS-I, the density profile exhibits a SF-Npp- 
Npp shell structure, while in the phase PS-II, the shell is 
in the form of SF-Npp. 

The analytical result of the density profile can be sum- 
marized as follows. In the region PS-I, we have 



n(r) 
n 




< x < x c 
x c < x < xq 
xq < x < X 



(51) 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 



r,=e B /e F 



FIG. 4: Global phase diagram for trapped 2D Fermi gas in 
the r\ — P plane. 




FIG. 5: The profiles for the total density n(r) (solid line) and 
magnetization m(r) (dashed line) in two cases, 77 = 0.2 and 
P = 0.5 in the region PS-I and 77 = 0.7 and P = 0.5 in the 
region PS-II. 
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and 



m(r) 
n 




26 

>.2 



< X < x c 

X c < X < Xq 
Xq < X < X 



(52) 



with no = cef, x = r/R^, Rt being the Thomas- Fermi 
radius of non-interacting gas defined in (|37p . and S = 
fife? being the real solution of the cubic equation S 3 — 
n 2 5/4:-Pn/2 = 0, 



5 = ^ 



Pr) 
1 



1/3 



(l+7) 1/3 + (l- 7 ) 



1/3 



(53) 



with 7 = y/1 -?y 4 /(108P 2 ). The scaled radii 
t c /Rt,Xo = tq/Rt and X = R/Rt are given by 



1 

2 25' 



Xq 

X 



V 



5, 



(54) 



A numerical sample for n = 0.2, P = 0.5 is shown in 
FigGUa). There is an interesting phenomenon which is 
different from that found in 3D: The magnetization pro- 
file m(r) exhibits a visible platform structure in the par- 
tially polarized normal shell in the region r c < r < r^. 
For partially polarized gas, the interaction may be im- 
portant jhke the finding around the unitary region in 
3D [HI [Hi l29j|. However, for the 2D system, since par- 
tially polarized normal shell appears only at small cou- 
pling where the effect of interaction is not important, our 
conclusion will not be qualitatively changed. We also ob- 
serve a density jump An at the critical radius r c . In the 
region PS-I, An is independent of the global polarization 
and depends only on the coupling strength, 



M 

An = nn = — e B . 

Ztt 



(55) 



Thus the experimental data for An can be used to extract 
the effective two-body binding energy cb- 
In the region PS-II, the density profile reads 



n(r) 



no 



2(1 -x 2 ) 



and 



m(r) 
n 



1 + (V2-1)VP- 



1 + (V2-1)VP-.t 2 



< x < x c 
x c < x < X 



(56) 



< x < x c 
x c < x < X 



The scaled radii x c and X now takes very simple form 



X r = VI -VP 



X 



1 + (V2-1)VP. 



(58) 



A numerical sample of the density profile for n — 0.7, P = 
0.5 is shown in Fig(5jb). It is very surprising that the 
density profile does not depend on the coupling param- 
eter n, but only on the global polarization P. As a re- 
sult, the ratio r c / R exhibits a universal behavior when 
77 > rji = 0.586, as shown in FigEl 




FIG. 6: The ratio of the superfluid radius to the cloud radius 
as a function of the global polarization P for n > ni = 0.586. 



Finally, two comments on our results should be made. 
The first is on the BCS-Leggett mean field theory. In this 
theory, the quantum fluctuation in the superfluid phase 
and the interaction in the partially polarized phase are 
totally neglected. The effect of interaction in the par- 
tially polarized phase may change the platform structure 
in the region r c < r < tq. However, since the three- 
shell structure appears in the weak coupling region, we 
expect this effect to be small. At very strong coupling, 
the correction in the superfluid phase due to quantum 
fluctuation should be important, and the universal be- 
havior in the region PS-II may be destroyed. Since the 
BEC region is reached at 77 > 2, we expect that our 
conclusion holds at the BCS side 771 < 77 < 2. The sec- 
ond comment is on the model we used. Recently, it is 
argued that the model we used is not sufficient to dis- 
cuss BCS-BEC crossover in quasi-2D Fermi gas due to 
the importance of dressed molecules [Hj]. However, from 
the study in [28|, this effect is important only at strong 
coupling (may be for 77 > 2). Obviously, the comparison 
of our prediction with the experimental data can tell us 
whether the quantum fluctuation, dressed molecules and 
other possible effects are important. 



VII. SUMMARY 

In summary, mean field phase structure of polarized 
Fermi gas in 2D is analytically investigated. In the nor- 
mal phase, the di-fermion bound states at high polar- 
ization are easier to survive than in the balanced Fermi 
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sea. In the BCS-Leggett mean field theory, the transi- 
tion from the unpolarized superfluid phase to the nor- 
mal phase is always of first order, and there exists no 
stable gapless superfluid phase. In the homogeneous 
system, we analytically determined the critical Zeeman 
fields and the critical population imbalance in the whole 
coupling parameter region. We found two critical Zee- 
man fields in the BCS-BEC crossover, and proved that 
the mixed superfluid-normal phase is the energetically fa- 
vored ground state. However, from recent Monte-Carlo 
simulations [TH [l(| , our mean field results may be only 
qualitatively correct in some parameter region, due to 
the importance of interactions in the normal phase. 

To compare our theoretical results with future exper- 
imental data, we have also calculated analytically the 
density profile for an imbalanced 2D Fermi gas confined 
in a harmonic trap. For balanced populations, the den- 



sity profiles for normal and superfluid matter are the 
same and can not be used as a signature of superflu- 
idity. For imbalanced populations, the density profile 
exhibits a shell structure, a superfluid core in the cen- 
ter and a normal shell outside. At small coupling, there 
exists a partially polarized normal shell and the density 
difference shows a platform structure. For large attrac- 
tion, however, the normal shell is fully polarized, and 
the density profile depends only on the global population 
imbalance. Our theoretical prediction can be examined 
in the future experiments on 2D ultracold Fermi gases, 
which can help us to see whether quantum fluctuations 
and other possible effects are important in determining 
the phase structure [Hj]. 
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